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Conventional approaches to thermal conductivity in itinerant systems neglect the contribution to 
thermal current due to interactions. We derive this contribution to the thermal current and show how 
it produces important corrections to the thermal conductivity in anisotropic superconductors. We 
discuss the possible relevance of these corrections for the interpretation of the thermal conductivity 
of anisotropic superconductors. 



Historically, understanding thermal conductivity has 
played an important role in the development of mate- 
rials theory. There has recently been a growth of in- 
terest in studying this basic property in strongly cor- 
related materials. 1,2 One of the recurrent consequences 
of strong interactions is the development of anisotropic 
gaps in the excitation spectrum: possible examples in- 
clude heavy fermion superconductors, Kondo insulators 
and the chevrel superconductor V3«Si. 3 ' 4 In such systems, 
the thermal conductivity has an important role to play 
in the elucidation of the gap symmetry. 

Pioneering work on the theory of thermal conductivity 
was carried out in the sixties 5 , by and large, the the- 
oretical approaches taken today are a direct application 
of this early work 8-10 . In this paper, we show how these 
classic approaches require modification to take account 
of new thermal conduction channels that are introduced 
by interactions. As part of this discussion, we shall show 
how thermal conductivity, like particle or charge conduc- 
tivity, can be regarded as a boundary condition response 
which is most directly computed as a response to a fic- 
titious gauge field. We illustrate this approach in the 
context of anisotropic superconductivity, showing how 
interaction contributions to the thermal current signif- 
icantly modify the flow of heat produced by excitations 
near a gap node. 

Thermal conductivity is directly related, via the Kubo 
formula, to thermal current fluctuations 
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There is an important distinction between heat and 
charge conductivity. Unlike charge density, the energy 
density does not commute with the interactions: this 
means that heat is directly transmitted by interactions, 
introducing new interaction contributions to the micro- 
scopic thermal current operator. The effect of these new 
thermal conduction channels is important when interac- 
tions severely modify the dispersion of the electrons, as 
in the case of a Mott insulator, or a superconductor with 
gap nodes. 11-13 

To derive the thermal current we appeal to Noether's 
theorem, which relates continuity of energy flow to the 
covariance of the action under co-ordinate transforma- 
tions in time t — > t' = t + cf>[x, t']. Consider an electronic 



system described by the Lagrangian 
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where H is the Hamiltonian and<9t= \{dt — dt) is the 
antisymmetrized time-derivative. By Noether's theorem, 
the change in the action S — J C(t)dt is the energy con- 
tinuity equation 



5S/S4>[x, t} = - [d t e(x, t) + V • j(x, t)] , 
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where e and j are the energy and thermal current den- 
sity respectively. Continuity of energy flow follows from 
the invariance of the action under co-ordinate transfor- 
mations 5S/8<p[x,t] = 0. Writing SS = J[SS/8<p]8<p(x,t) 
and integrating by parts, we find 
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dtdx 
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so that e(x) — 8S/5<j), while the heat current is j = 
SS/SVcj). The calculation of the thermal current requires 
special care because interactions are non-local and all 
higher derivatives of the fields must be taken into ac- 
count when taking the functional derivative. This im- 
portant point was overlooked in the classic treatment by 
Langer 6 . We can take these effects into account by noting 
that the derivative operators inside the action are evalu- 
ated at constant time t, so that under the transformation 
t — ► t', they acquire a covariant form 
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where we denote Co = i dt> ■ The field A = V</> adjusts 
for the fact that spatial derivative V = V|f is taken at 
constant t' . It follows that the energy current is 



j = 8S/SV</) = -dH[A]/dA, 



(6) 



where H[A] is the "gauged" Hamiltonian with the re- 
placement k — > k — Auj in both the kinetic and interac- 
tion terms. We see that A plays the role of a "fictitious" 
gauge field conjugate to the thermal current. 
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To illustrate this procedure, consider a fluid of elec- 
trons with kinetic energy eg and an exchange interaction 
Jg, where 

H[A] = ]T [^(k-A^k + \°-k ■ J(k-A^k\ > (7) 



where Bg is the Fourier transform of the spin density at 
wavevector k . Differentiating with respect to A gives 
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The second term reflects the additional heat flow created 
by the exchange interaction. In an antiferromagnetic 
Mott insulator, it is this term which is responsible for the 
conduction of heat by spin- waves. The approach illus- 
trated here can be used on any Hamiltonian, without the 
need to develop the equations of motion for that partic- 
ular model. Past efforts to compute the thermal current 
contribution from interactions have adopted an equation 
of motion approach on a case- by-case basis. This yields 
long expressions for the thermal current, where the time 
derivatives are expanded in their full glory. 11 ' 12 

The usefulness of the gauge-theoretic derivation of the 
thermal current lies partly in its invariancc under the 
renormalization group. This means that it can be di- 
rectly applied to the effective Lagrangian that describes 
the low energy physics of an interacting system. In any 
case where the low energy physics is described by weakly 
interacting quasiparticles with energy Eg, application of 
the same procedure yields 
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where ng 



alar is the quasiparticle number operator 

k K 

and we have used the equation of motion to make the 
last substitution. When the mass renormalization of the 
electrons is highly anisotropic, the quasiparticle thermal 
current contains a large interaction component. Micro- 
scopic calculations that ignore these contributions fail to 
recover the correct quasiparticle description of the heat 
current at low temperatures. 

A particularly important illustration of this effect oc- 
curs in an anisotropic superconductor. In this case, the 
BCS Hamiltonian takes the form 

ff[I] = l/2^^^ 
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where Ag is the anisotropic gap function, \I/t is the four- 
component Nambu spinor for the electrons. The thermal 
current operator is then 





FIG. 1. Constant energy contours surrounding a line node 
in a "polar" superconductor. The thermal current operator is 
normal to the energy contours and points radially outwards 
from the gap node. The dense line indicates position of Fermi 
surface and the dark region represents the line node, which 
lies perpendicular to the plane of the paper. 
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The off-diagonal terms in this expression derive from the 
interaction contribution to the thermal current. In the 
presence of weak scattering, we can use the equation of 
motion to replace 

v fe ^^ - \{e % , v s e s ] = %v g %i, (12) 

where Eg = {el + A?-) 1 / 2 . This recovers the quasipar- 
ticle form (9). Note however, that the Fermi veloc- 
ity is replaced by the the quasiparticle group velocity 
Vg = WgEg, which has important components parallel 
to the Fermi surface. These can significantly affect the 
anisotropy of the thermal conductivity. We now use (11) 
and apply the Kubo formula, to obtain 

^ = ^E /~ dujuj2 (-£) A ^' w )' ( 13 ) 



where 
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Ag(w) is the matrix spectral function and / is the Fermi 
function. Carrying out the trace we find 
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where T(u) is the quasiparticle scattering rate and 
N{w){VV^ = Y,^ j ^-E % ) 



(16) 



describes the quasiparticle velocity distribution. 

To show how the interaction contribution to the ther- 
mal current modifies the thermal conductivity, consider 
the case of a "polar" superconductor, with a spherical 
Fermi surface and gap function |A^| 2 = Aflfcf, which 
leads to a line of gap nodes in the basal plane k z = 0. 
The kinetic contribution to the thermal current operator 
always lies normal to the electron Fermi surface, and if we 
omit the effects of interaction, the low temperature ther- 
mal conductivity lies predominantly in the basal plane. 
In the correct calculation, the current carried by each ex- 
citation lies in the direction of the quasiparticle velocity 
, which lies normal to the surfaces of constant energy 
(Fig. 1). The line- node in the basal plane is surrounded 
by elliptical surfaces of constant energy and consequently, 
in the vicinity of the gap node there are quasiparticles 
which propagate perpendicular to the basal plane, giving 
rise to a increased thermal conductivity in this direction. 
At low temperatures the thermal conductivity in both 
directions is proportional to T 2 , with a finite ratio at 
T = 



1 (Ao 
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Without interaction effects, this ratio vanishes at T = 0. 

In Figure 2 we show how the anisotropy in the ther- 
mal conductivity varies as the ratio between the maxi- 
mum gap and the Fermi energy is increased from zero. 
Two effects of interest are observed: (i) the anisotropy 
is reduced by these interaction effects; (ii) the ther- 
mal conductivity is enhanced in the region close to the 
transition temperature. These effects become important 
when the transition temperature is a substantial frac- 
tion of the Fermi energy. Such considerations may be 
particularly relevant to heavy fermion superconductors, 
of which two cases are particularly worthy of mention: 
UBeis and UPt%. In both cases, NMR measurements 
suggest the presence of line nodes. In the former case, 
superconductivity develops before a Fermi liquid regime 
has been established, so it is very likely that A/ep is 
large. In the latter case, the conventional wisdom is that 
A/ep ~ 0.05 — 0.1. However, measurements of the ther- 
mal conductivity on very pure samples of this material 1 
show that k zz I k xx ~ constant at low temperatures; they 
also display a marked peak in the ratio k sc /k u near T c . 
Existing approaches 2 attribute these effects to a com- 
bination of both elastic and inelastic scattering. While 
such an explanation is surely feasible, an alternative ex- 
planation might be obtained by assuming a larger ratio 
A/ep, and attributing these features to the interaction 
contribution to the thermal current that has been hith- 
erto ignored. 
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FIG. 2. Ratio between the thermal conductivity perpen- 
dicular and parallel to the basal plane as a function of 
the normalized temperature for different values of the ratio 
r = A /e F . 

The main point of this paper has been to emphasize 
that interactions actually modify the thermal current, 
which can have important effects on the thermal conduc- 
tivity. As an illustration of these ideas, we have consid- 
ered the effect of gap anisotropy in polar superconduc- 
tors. Our basic approach can be applied to other gap 
anisotropic systems. One very interesting case in this 
respect may be the small gap Kondo insulators, where 
the insulating gap appears to exhibit a point node 14 . 
Similar considerations can also be applied to the mo- 
mentum current, a point which is important in the in- 
terpretation of measurements of the ultrasound attenu- 
ation in strongly correlated systems. This could explain 
why early ultrasound measurements in UBeis 15 found 
no anisotropy, despite the clear suggestion of line nodes 
from NMR measurements. 16 

A secondary aspect of our work concerns the observa- 
tion that thermal currents can be treated using a gauge 
theoretic approach. Just as momentum currents, or 
stress (a) in a material result from the gauge field we 
call strain (u), a = Ru, energy currents result from a 
"temporal strain" which we represent by the gauge field 
A. If we calculate the heat current response to this ficti- 
tious field 



j T (u) = -Q(u)A(u>), 



(18) 



we find that the "London Kernel" is directly related 
to the thermal conductivity by the relation Q(u>) = 
1(jjTk,{ui), from which we deduce that 
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In other words, A is the gauge field conjugate to thermal 
gradients. By analogy with the case of momentum cur- 
rents, the quantity A may be regarded as a sheer in time. 
Just as broken spatial translation symmetry enables a 
material to support a persistent stress or "momentum 
superflow" , we are led to speculate that if it were possi- 
ble to produce a state of condensed matter with sponta- 
neously broken time translation symmetry, such a system 
would exhibit "heat superflow" , manifested by an infinite 
thermal conductivity. 
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